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Abstract 

This note provides a simple result showing, under suitable technical assumptions, that if a system 
E adapts to a class of external signals U, in the sense of regulation against disturbances or tracking 
signals in U, then E must necessarily contain a subsystem which is capable of generating all the 
signals in U. It is not assumed that regulation is robust, nor is there a prior requirement for the 
system to be partitioned into separate plant and controller components. Instead, one assumes that 
a "signal detection" property holds. 
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1 Introduction 

Suppose that one knows that a certain system E regulates against all those external input signals u 
which belong to a predetermined class U of time- functions. (Input signals u are often thought of as 
disturbances to be rejected or signals to be tracked, depending on the application.) In biology, one 
often uses the term adaptation for this property. It means that a certain quantity y(t) associated to the 
system, called its output (also called a regulated variable or an error) has the property that y(t) — > 
as t — > oo whenever the system is subject to an input signal from the class U (Figure . The internal 



u(-) eW 



y(t) -► 



Figure 1: Given System, Regulated Output y(t) when Inputs in U 



model principle (IMP) states, roughly, that the system E necessarily must contain a subsystem Etm 
which can itself generate all disturbances in the class U. The terminology arises when thinking of Etm 
as a "model" of a system which generates the external signals. 

For example, if y(t) — > as t — > oo whenever the system is subject to any external constant signal 
(i.e., the class U consists of all constant functions), then the system E must contain a subsystem Etm 
which generates all constant signals (typically an integrator, since constant signals are generated by the 
differential equation ii = 0). Of course, the choice of y = as the "adaptation value" is merely a matter 
of convention; by means of a change of variables, one may always reduce a given regulation objective 
"y(t) — > yo" where yo is some predetermined value, to the special case yo = 0. 

*E-mail: sontag@control.rutgers.edu. Supported in part by US Air Force Grant F49620-01-1-0063 and NIH Grant 
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In addition, the IMP specifies that, in an appropriate sense, the subsystem Ejm must only have y 
as its external input, receiving no other direct information from other parts of the system nor the input 
signal u. One intuitive interpretation is that Sim generates its "best guess" of the external input u based 
on how far the output y is from zero. Pictorially, if we have the situation shown in Figure ^ then there 
must be a decomposition of the system E into two parts, as shown in Figure [2 where the system Sim 



u(-) eU 




y(t) -> 



Figure 2: Decomposition of S into E and Sim, the Latter Driven by y(t) 



(with y = 0) is capable of reproducing all the functions in 1A. 

The internal model principle originates in the biological cybernetics literature. But, as with any 
"principle" in control theory (like dynamic programming, the maximum principle, etc.) and more 
generally in mathematics, the IMP is not a theorem but rather a "mold" for many possible theorems, 
each of which will hold under appropriate technical assumptions, and whose conclusions will depend 
upon the precise meaning of "class of external signals" , "reproducing all functions" , and so on. 

The best known instance of an internal model theorem is due to Francis and Wonham, who in a series 
of beautiful and deep papers in the mid 1970s proved a theorem for linear systems which showed, in 
essence, that structurally stable or "robust" adaptation forces the existence of embedded internal models. 
Partial generalizations of their work to nonlinear systems were later obtained by Wonham and Hepburn, 
see n,5Hl,CHjj|2]-|Zj- The Francis/ Wonham theory applies to systems S which are already partitioned 
into a "plant" plus a "controller" . The robustness assumption amounts to the requirement that the 
given controller should perform appropriately (in the sense that the regulation objective y(t) — > 
is achieved) even when the plant subsystem - but most definitely not the controller subsystem - is 
arbitrarily perturbed. The conclusion is that the controller is driven by y and incorporates a model of 
the external signals. That some additional condition - such as structural stability - must be imposed 
is obvious, since the system E which simply outputs y = for every possible input signal u does not 
contain any subsystem generating the signals u. We will impose instead a condition which amounts to 
a signal detection property: the output must reflect sudden changes in the input. 

Note that this type of objective is very different from what would be typical in control design: in 
the latter field, one would ideally not even notice disturbances (for instance a change in the road grade, 
in an automobile's cruise-control system, or a bump in the road, in an active suspension system) In 
contrast, in biological applications, signal detection is often an objective, to be followed by a return to 
default values. This subtle difference in desired behaviors, while dealing with what are otherwise similar 
problems, is characteristic of many applications of control-theory ideas in biology. 

Recent work in molecular biology, cf. J7|, nas suggested that the IMP could help guide experimen- 
talists and modelers: if certain characteristics of a system adapt to signals in a given class (in all the 
examples so far, constant inputs, such as for instance y(t) = the relative "activity" of enzymes con- 
trolling motors in E.coli chemotaxis, with respect to u(t) = concentration of extracellular ligand, but 
similar considerations may apply to periodic inputs and circadian clocks as internal models of day /night 
periods) then the IMP could, in principle, help distinguish among mathematical models which do or do 
not contain internal models. 

With a view toward such biological applications, it is desirable to have available a theorem which 
(a) applies to nonlinear systems E, at least under reasonable technical assumptions, and (b) does not 
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require the system E to be split between "plant" and "controller" subsystems, nor (c) requires structural 
stability (robustness) in the sense of the Francis/ Wonham theory (which would imply, in the case of the 
E.coli motor control network, that the system should perfectly adapt even if there are arbitrary direct 
connections between the external ligands and the motor signals, a matter which seems difficult to check 
experimentally), and relies instead upon a signal detection property. We present one very elementary 
and self-contained such result in this note. It basically just picks from and "repackages" some of the 
basic concepts and techniques developed by previous researchers for the same problems, in particular: 
the use of differential geometric techniques and "output-zeroing" sets (^E|, EHD; dynamical 
systems notions like omega-limit sets (0, 0,^01)) an d system decompositions motivated by the Center 
Manifold Theorem ([5], 0,0, EU)- Isidori's excellent textbooks ^2;E2| should be consulted for a far 
deeper discussion of many of the issues raised here. 

Precise mathematical definitions are provided in Section On the other hand, since the linear 
version of the result is very easy to explain, we sketch that case first. (The discussion assumes some 
familiarity with frequency domain techniques, and may be skipped without loss of continuity.) 

1.1 Linear Case 

Let us denote by S the transfer function of the system E: if y is the output produced when E starts at 
the zero initial state and is fed input u, then the relation y(s) — S(s)u(s) holds between the Laplace 
transforms y(s), u(s) of the output and input. One expresses S(s) = |W as the quotient of two 
relatively prime polynomials, with the degree of p less than the degree of q. (An equivalent discussion 
using differential operators instead of Laplace transforms is also possible, see e.g. Section 6.7 in |14|.) The 
first observation, a well-known fact in systems theory, is that the zeroes of p can be viewed, alternatively, 
as poles of a feedback subsystem. To see this, we assume that p is not identically zero, and divide the 
polynomial q by p, obtaining q = ap + &, where b is some polynomial of degree less than p. Now, as 
the algebraic equality y = -u is equivalent to y = -(u — -y), we conclude that the system E can be 
decomposed as in Figure^ For example, if s — is a zero of S (that is, is a root of p), which amounts 
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Figure 3: System Equivalent to E: Closed-loop Zeros are Feedback Poles 

to the property that constant signals get differentiated by E (the "DC gain" of E is zero), then the factor 
1/s appears in the feedback box b/p, and can be interpreted as an integrator of the output y. 

We will show that, in general, the subsystem with transfer function ^ models all inputs which E 
adapts to. Let us suppose that the class U of inputs can be described as the set of all possible solutions 
of a fixed linear differential equation 

u (( \t)+b 1 u (t - 1) {t) + ... + b i - 1 u'(t) + b e u(t) = 

for some integer £, and which has no stable modes. (Stable modes, giving components of u which converge 
to zero, are less interesting, since they do not represent persistent disturbances.) We view these signals u 
as the outputs of an "exosystem" T which is obtained by rewriting the differential equation as a system 
of £ first order equations. Figure^Jshows a cascade consisting of the original system E and the exosystem 
F which generates the inputs in U. (If, for example, U — constant inputs, then one would let T be the 
system with equation w = and output u — w, and for each initial condition w(0) — u°, one obtains 
a different constant output u(t) = u°.) The regulation objective is now simply that y(t) — > for all 
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Figure 4: Exosystem and System in Cascade 



possible initial conditions of the composite system, i.e. for all initial conditions of the original system S 
and all initial conditions of the exosystem T, the latter corresponding to all possible inputs in U being 
fed to E. 

Next, we reformulate this regulation property by adding an external input v(-) to the exosystem, and 
requiring now that y(t) — ► for all possible stable inputs (v(t) — > as t — > oo) but only when starting 
from the zero initial state. (Such replacements of initial states by stable forcing inputs - assuming 
natural controllability/observability conditions - are elementary exercises in linear systems theory, see 
e.g. the proof of Theorem 33 in |14j . ) In other words, we have now the situation illustrated by FigureEl 
We denote by G the transfer function of the exosystem T: G = —, where 



,(t) - 




Figure 5: Exosystem and System, Forced by Stable Inputs 



tt(s) = s i + b l s i 1 + . . . + b ( -is + b e . 

To see that the subsystem with transfer function b/p includes an internal model of T, we argue as follows. 
The regulation property for the cascade in Figure|S]means that the product rational function GS is stable 
(all poles have negative real parts), while the assumption that G had no stable modes means that all the 
poles of G (i.e, the roots of the polynomial it) have nonnegative real parts. Therefore, these poles must 
be canceled in the product GS] in other words, S must have among its zeroes all the poles of G, so that 
we can write p = npo for some polynomial pq. Thus b/p — b/(Ttpo). One may now factor b = &1&2 in such 
a way that the degree of &2 is less than the degree of %, so that b/p = (fri/poXfe/V) and now the system 
with transfer function b/p can be written itself in the cascade form in Figure El The subsystem with 
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Figure 6: Decomposition of - 



transfer function 62/^ generates all the inputs in U, since one may write a set of differential equations 
for it which is exactly the same as for the exosystem T, changing only the output mapping ("controller 
form" realization). 

Since the tools of transfer functions are not available for nonlinear systems, a different approach is 
required in general. 
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2 Definitions and Statement of Result 



We consider single-input single-output systems S, affine in inputs: 

x(t) = f(x(t)) + u(t)g(x(t)) , y(t) = h(x(t)) (1) 

(dot indicates derivative with respect to time, and the arguments t will be omitted from now on; see |14| 
for general definitions and properties of systems with inputs). Here x(t), u(t), and y(t) represent the 
state, input, and output at time t, f and g are smooth vector fields on R™ (n is the dimension of the 
system), h is a scalar smooth function R n — > R, and /(0) = h(Q) = 0. (Several assumptions on / and g 
will be made later.) A special case is that of linear systems 

x = Ax + ub , y = cx (2) 

where A is an n x n matrix, b is a column n-vector, and c is a row n-vector. 

Suppose given a class U of functions [0, oo) — > R (such as for example the set of all constant functions). 
We say that E adapts to inputs in U (a more appropriate technical control-theoretic term would be 
"asymptotically rejects disturbances in W) if the following property holds: for each u £ U and each 
initial state a; £ R n , the solution of Q with initial condition x(0) = x° exists for all t > and is 
bounded, and the corresponding output y(t) = h(x(t)) converges to zero as t — > oo. 

We will say that S contains an output-driven internal model oflA if there is a change of coordinates 
which brings the equations Q into the following block form: 

Zl = fi(zi,z 2 ) +ug 1 (z 1 ,z 2 ) 

h = h(y,z2) (3) 

V = 

(the subsystems with variables z\ and z 2 correspond respectively to So and Sim m Figure^, and in 
addition the subsystem with state variables z 2 is capable of generating all functions in IA, meaning the 
following property: there is some scalar function ^(22) so that, for each possible u £ U, there is some 
solution of 

i 2 = / 2 (0,z 2 ) (4) 

which satisfies <p(z2 (t)) = u(t). 

The precise meaning of "change of coordinates" is as follows. There must exist an integer r < n, 
diffcrentiable manifolds Z\ and Z 2 of dimensions r and n — r respectively, a smooth function k : Z\ — > R, 
vector fields F and G on Z\ x Z 2 which take the partitioned form 

F= ( fi(zi,z 2 ) \ q _ fgi{zi,z 2 )\ 
\.f 2 (Mzi),z 2 )J ' \ J 

and a diffeomorphism $ : R" — > Z% x Z 2 , such that 

$*(x)f(x) = F(<t>(x)), *.{x)g(x) =G($(x)), k($i(x)) = h(x) 

for all x £ R™, where $i is the Zi-component of $ and star indicates Jacobian. 

Our result will hold under additional conditions on the vector fields defining the system. The first 
condition is the fundamental one from an intuitive point of view, namely that the system is able to 
detect changes in the input signal: 



Assumption 1: a uniform relative degree exists. 



This means that there exists some positive integer r such that 

LgLjh = Vfc<r-1 
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and 

LgL^hix) ^ Vie R n 

where, as usual, Lxh indicates the directional derivative ("Lie derivative") of a function h along the 
direction of the vector field X, that is (Lxh)(x) = Vh(x) ■ X(x). The integer r, if it exists, is called the 
relative degree of E. It is possible to prove (see [llj l that when r exists, necessarily r < n. 

For a linear system (J2J, existence of a relative degree amounts to simply asking that cA % b is nonzero 
for some i, or equivalently that the transfer function c{sl — A)~ 1 b is not identically zero. For general 
systems (JTJ, the assumption is equivalent to the statement that the output derivatives y^(t) must be 
independent of the value of the input at time t, for all k < r, but that y^ r '(t) = b(x(t)) + a(x(t))u(t) 
for some function a(x) which is everywhere nonzero (so that the system can be "inverted" to obtain the 
instantaneous value u(t) from instantaneous derivatives). See also for a discussion of the character- 
ization of r in terms of smoothness of outputs when inputs are discontinuous (change detection) . 

The next two conditions are of a technical nature. They are automatically satisfied for linear systems. 
For nonlinear systems, we need such conditions in order to guarantee the existence of a change of variables 
exhibiting the system Sim- (See Remark 13.21 for ways of weakening these assumptions.) We are guided 
by conditions which appear in Isidori's book 

Assuming that the degree is r, we introduce the following vector fields: 

g(x) = L L r-i h ^ 9( x ) i f( x ) = f(x) ~ ^L r f h(x)Jg(x) , n := ad^g, i = 1, . . . r , 

where adx is the operator adxY = [X, Y] = Lie bracket of the vector fields X and Y. Recall that a 
vector field X is said to be complete if the solution of the initial value problem x — X(x), x(0) — x° is 
defined for all tgl, for any initial state x , and that two vector fields X and Y are said to commute if 
[X, Y] = 0. The assumptions are: 



Assumption 2: Tj is complete, for i = 1, . . . , r. 



Assumption 3: the vector fields r, commute with each other. 



Finally, we must define the allowed classes of inputs U. As usual in control theory (see also the 
discussion in Section ll.lfl . we will assume that inputs are generated by exosystems. That is, there is 
given a system T: 

w = Q(w), u = 6(w) (5) 

(let us say evolving on some differentiable manifold, Q a smooth vector field, and 9 a real-valued smooth 
function, although far less than smoothness is needed) such that the input class IA consists exactly 
of those inputs u(t) = Q(w(t)), t > 0, for all possible solutions of w = Q(w). For example, if we are 
interested in constant signals, we pick w = 0, u = w and if we are interested in sinusoidals with frequency 
ui then we use xi = x%, x% — —ui 2 xi, u = x\. It is by now standard in nonlinear studies of necessary 
conditions for regulation to impose conditions on omega limits sets for trajectories of the exosystem, 
see 0,0; we will follow the approach in 50j-C3 an d assume that the exosystem is Poisson-stable: for 
every state w°, the solution w(-) of w = Q(w), w(0) — w° is defined for all t > and it satisfies that 
w° is in the omega-limit set of w, that is, there is a sequence of times i; — > oo such that the sequence 
w(ti) converges to w° as t — > oo. This means that the exosystem is almost-periodic in the sense that 
trajectories keep returning to neighborhoods of the initial state. 

This theorem is proved in Section 

Theorem 1 If Assumptions 1-3 hold and the system S adapts to inputs in a class U generated by a 
Poisson-stable exosystem, then S contains an output-driven internal model ofU. 
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2.1 An Example 



As an example, consider the model for E.coli chemotaxis adaptation to constant inputs given in 
Section 2.2. Letting x\ = R and X2 = RL be the concentrations of unbound and bound receptors re- 
spectively, and taking the external ligand concentration u = L as input, we have the following equations: 

x\ = a% — 02X1 + 0,3X2 — CL4X1U 

for suitable positive constants a%, . . . , a§. In terms of vector fields, 

, _ fax — a 2 xi + a3X 2 \ _ f— (14XX 

\ a 5 - a 6 x 2 J ' 9 \ a^x-i 

and, still as in HJ, we take as output y the difference between the total concentration of active receptors 
and a steady state level of this activity. In terms of the notations used here, and up to multiplication 
by a suitable constant, this amounts to the following choice: 

h(x) = A Q - A = [ax + a 5 ] - [a 2 x x + (a 6 - a 3 )x 2 ] . 

We note that L g h = Dxx, where D = a 2 ax + (03 — ae)ax- Except in the accidental case when this 
constant D vanishes (in terms of the notations in 9 , D — k-ilTk r {ax — ct 2 ), so D can only vanish if 
ai = a 2 ), we have that L g h(x) ^ for all x {x\ > 0, as it represents a concentration), and so it follows 
that E has well-defined relative degree r = 1. Moreover, n = g is a constant vector, so Assumptions 2 
and 3 hold as well. 

A minor technicality concerns the assumptions that our systems (QJ evolve in all of Euclidean state 
space (not just xi > 0) and that /(0) = h(0) = 0. However, this is just a matter of picking the right 
coordinates. Notice that / vanishes at x° = {x\,x 2 ), where x\ = (axa& + a3a5)/(a 2 ae) and x 2 = a^/ae, 
and h vanishes at x° too. In order to fit into the general theory, one simply changes variables, mapping 
the positive orthant into all of M 2 and x° into the origin by means of x\ — lnxj — lnx". (Of course, 
there is no need to actually perform the coordinate change, since conditions expressed in terms of Lie 
derivatives are covariant.) 

Finally, letting B := x\ + x 2 (as done in one obtains a system of equations in terms of the new 
variables A and B, for which B = y. This last equation represents an integrator (internal model of a 
system which produces constant inputs) driven by the output y. (Of course, there is no point in applying 
the theorem, since once that the model is given we can find the internal model explicitly.) 



3 Proof of Theorem [T] 

Suppose that the system E adapts to inputs in U, which are produced by a Poisson stable exosystem 
r. We consider the interconnected system consisting of the cascade of F and E, as shown in Figure 0J 
namely: 

w = Q{w) , , 

x = f{x) + e(w)g(x) {l} 

and let Z denote the set consisting of those states x of E for which h(x) = (the "output-zeroing" 
subset). 

Lemma 3.1 For each w° there is some solution a — (w(-), #(•)) of the composite system {7J such that 
w(0) = w° and x(t) € Z for all t > 0. 

Proof. We start by picking an arbitrary solution cr = (w(-),x(-)) of the composite system (J7J such that 
u>(0) = w°, and let £1 = f2 + [cro] be the omega-limit set of this trajectory. We claim that, for each point 
(w,£) G (we partition coordinates into those for T and E) it must be the case that £ € Z. Indeed, by 
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definition of SI there is some sequence of times tj — > oo such that a;(i,) — » £. Since h(x(U)) — > because 
of the adaptation property and /i is continuous, it follows that h(£) = 0, as claimed. Next, we claim that 
there is some a; such that (w°, x°) £ ft. To see this, we first pick a sequence of times U — > oo such that 
w(i,) — » w° (Poisson stability is used here); as {x(ti)} is bounded, we may pick a subsequence U j of the 
ti so that x(ttj) — ► a: for some and this proves that (w°,x°) S 0, as wanted. 



Finally, we let a be the solution a = (w(-),x(-)) of the composite system (JJJ) for which w(0) 



w 







and a:(0) = a; , where x° is so that (w°,x°) G SI. Omega-limit sets are invariant, so a(t) G SI for all 
t > 0, and we already proved that this last property implies that x(t) £Z. I 

Proposition 9.1.1 in 11, shows that there is a global diffeomorphism $ so that, in the new coordinates, 
the system S takes the form shown in Display 0- Moreover, the subsystem described by z\ evolves in 
W and, using coordinates z\ = (Ci, ■ ■ • , Cr)j the equations for z\ can be written as follows: 



Ci = c 2 

Cr— 1 Cr 

Cr = b(z 1 ,z 2 ) + a(z 1 ,z 2 )u 



(8) 



where the output is y — k(zi) = Ci and a, b are smooth functions with a(z) — L g L r ^ 1 h(<& 1 (z)) ^ for 
all z. Wc let 

a{0, z 2 ) 

and show that for each possible u GU there is some solution of (@J which satisfies ip(z2(t)) = u(t). 

We pick to such that u(t) = 9(w(t)) and w(0) = w°, and view the interconnection Q of T and S in 
terms of the coordinate change given by $ on S: 

w = Q{w) 

zi = fi(zi,z 2 ) + 9(w)gi(z 1 ,z 2 ) 
Z2 = f2{y,z 2 ). 

Lemma 13.11 gives us the existence of a solution a — (w(-), zi(-), z 2 (-)) such that Q(w(t)) = u(t) and 
(i(t) = 0. Because of the form (JSJ of the zi-subsystem, this implies that zi(t) = and that ( r {t) = 0. 
Thus, along the solution a one has 6(0, z 2 (t)) + a(0, z 2 (t))u(t) = 0, and this is precisely what we wished 
to prove. I 

Remark 3.2 Assumptions 2 and 3 are automatically satisfied for linear systems, since the vector fields 
Ti are all constant, so that they are indeed complete and pairwise commutative. For general nonlinear 
systems, these assumptions, especially 3, are quite strong. Weaker conditions may be given, if one is 
merely interested in a local result, or if one is willing to accept a subsystem Sim which is driven by not 
just y but also several derivatives of y. Indeed, assuming merely a well-defined relative degree around 
a given point a; , we obtain a decomposition as in © (see Section 4.3 in except that the change 
of coordinates is now only valid in a neighborhood of a; , and f 2 now depends on (zi, z 2 ) (as opposed 
to (y, z 2 ). Note that, from the form JHJ, Z\ is the vector consisting of the derivatives y,y', ■ ■ - y^- If 
condition 2 holds, but 3 does not, then a global decomposition is possible, but f 2 still depends on 
derivatives of y (cf. Proposition 9.1.1). □ 



A Remark on Subsystems 



We expressed our theorem in terms of the existence of solutions which reproduce all inputs. Under 
additional and stronger hypotheses, one could also obtain an actual embedding of the exosystem in the 
internal model Eim ■ A full nonlinear version would involve abstract quotients of systems under suitable 
equivalence relations, and may follow along the lines of the work in 0] (based on 0). However, the 
necessary steps are easy to understand and prove in the case of linear systems. We start by showing the 
following elementary fact from linear systems theory 
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Lemma 3.3 Suppose given an observable linear system w = Qw, y = 9w and another linear system 
Z2 = Fz2 + Gy, u — ipz, and assume that for each w° there is some z° such that ipe tF z° = 6e t( ^w for 
all t > 0. Then, the matrix F is similar to a matrix with this block structure: 



(9) 



Proof. We first assume that the pair (F, tp) is observable, and claim that for each uP there is a unique z° 
such that ipe tF z° = 8e t( ^w°. This is because ipe tF z° = ipe tF z 1 implies z° = z 1 (observability). So we can 
define a map T : w° i— > z°. This map is one-to-one, by observability of the pair (Q,8). It is also linear, 
since 9e tQ (aw° + w 1 ) = a6e tQ w° + de^w 1 = a(pe tF Tw° + ipe^Tw 1 = cpe tF (aTw° + Tw 1 ) means that 
aw° + w 1 t— > aTw° + Tw 1 . It also satisfies FT = TQ, since taking derivatives in ipe tF Tw° = 9e t( ^w° 
gives (pe tF FTw° = 9e t( ^Qw° which means that Qw° i— > FTw°. Thus, on some invariant subspace (the 
range of T), F can be written as Q, which means that we can write F up to similarity in the form 
Q ^ 
*, 



Since F is similar to its transpose, and Q is similar to its transpose, F is also similar to 



a matrix in the form ( J. An observability decomposition (JUj Chapter 6) then reduces to the 



observable case. 



Without loss of generality, one may assume that linear exosystems are observable (there always 
exists an observable equivalent). We now apply Lemma i3.3l to the exosystem and the internal model 
Sim, assumed linear. There results a change of variables for Sjm so that, in the new variables, a subset 
£ of the variables z 2 of Sim, corresponding to the first block in ©, evolves according to an equation of 
the form £ = QC, + by, for a suitable vector b. This provides the desired embedding of the exosystem in 
the internal model. 
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